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Abstract. The key condition A3w of Ma, Trudinger and Wang for regularity 
of optimal transportation maps is implied by the nonnegativity of a pseudo- 
Riemannian curvature — which we call cross-curvature — induced by the 
transportation cost. For the Riemannian distance squared cost, it is shown that 

(1) cross-curvature nonnegativity is preserved for products of two manifolds; 

(2) both A3w and cross-curvature nonnegativity are inherited by Riemannian 
submersions; and (3) the n-dimensional round sphere satisfies cross-curvature 
nonnegativity. From these results, a large new class of Riemannian manifolds 
satisfying cross-curvature nonnegativity (thus A3w) is obtained, including 
many whose sectional curvature is far from constant. All known obstructions 
to the regularity of optimal maps are absent from these manifolds, making 
them a class for which it is natural to conjecture that regularity holds. This 
conjecture is confirmed for complex projective space CP". 



This paper addresses questions in optimal transportation theory and Riemannian 
geometry. For a general introduction to these subjects we refer to the books by 
Villani [Vlj |V2j for optimal transport theory and the book by Cheeger and Ebin 
CE for Riemannian geometry. 

1.1. Background: optimal transport and pseudo-Riemannian geometry. 

In optimal transportation theory one is interested in phenomena which occur when 
moving mass distributions so as to minimize the transportation cost. Mathemat- 
ically, there are source and target domains, M, M, two differential manifolds 
equipped with a lower semi-continuous cost function c : M x M — ► R U {oo}. 
Given two positive Borel probability measures p, p on M, M, respectively, one 
wants to understand the optimal map F : M — > M , which minimizes the average 
cost 
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among all such maps pushing- forward p to p, denoted F#p = p and meaning that 
p(F- l {E))=p{E) V Borel E C M. 

Particular attention has been devoted to the case M = M, a complete Riemann- 
ian manifold with the cost c = ^dist 2 , where dist denotes the Riemannian distance 
function. Existence and uniqueness of optimal maps in this case is well known 
due to the work of Brenier [Br] and McCann [M] (and Figalli [F] for noncompact 
manifolds), under the condition that p doesn't charge lower dimensional submani- 
folds (see also |CNMj ). Under suitable hypotheses, regularity (C° jC a fC°°) of such 
optimal maps is known for Euclidean space M = R" by results of Delanoe |D1| . 
Caffarelli |Cal| |Ca2] , and Urbas [U] , and for fiat [Coj and near flat |D2j manifolds 
by Cordero-Erausquin and Delanoe. Beyond the flat case, Loeper |L2j deduced reg- 
ularity on the round sphere, by combining his own breakthroughs with pioneering 
results of Ma, Trudinger and Wang [MTWj [TWlj |TW2| concerning regularity of 
optimal maps for general cost functions. Loeper's result is simplified in the work 
of present authors KM where we give an elementary and direct proof of a crucial 
maximum principle deduced by Loeper from Trudinger and Wang's theory. We re- 
ferred to Loeper's maximum principle with the acronym DASM (see Theorem l2.7|) . 
and extended it to the manifold case. In the course of deriving new results, our 
method was employed and further developed by Figalli and Villani |FV] . Figalli 
and Rifford [FRj, Loeper and Villani [LV], and Villani [V2] [V3] . 

Our contributions in [KM] are based on a pseudo-Riemannian geometric struc- 
ture h on M x M induced by the transportation cost. Namely when dim M = dim M 
this pseudo- metric h is defined on N C M x M with c £ C A (N) as the following 
symmetric bilinear form on TM TM: 

where D and D denote the differentials along M and M, respectively. For non- 
degeneracy of h we assume that DDc and its adjoint DDc are non-degenerate 
(Ma, Trudinger and Wang's condition A2: see Section [2]). This is automatically 
true for the Riemannian distance squared cost, away from the cut-locus, since 
DDc is a matrix of independent Jacobi fields. This pseudo-metric h geometrizes 
the regularity theory of optimal maps, by recasting the key cost hypothesis A3w 
of Ma, Trudinger and Wang [MTWj [TWlj for regularity as the non-negativity of 
certain pseudo-Riemannian sectional curvatures of h. To explain this condition 
more precisely, we need 

Definition 1.1 (cross-curvature). Let (x, x) S N C M x M. For each p(&p € 
T(x,x)N = T X M © T S M, the cross-curvature of p and p is defined as 

(1.2) ctossJ^Vp) =Rh((p®0) A (0©p),(j)fflO)A (0 9p)) 

where Rh is the Riemann curvature operator of the pseudo-metric h. We drop the 
superscript (N, h) and the subscript (ir, x) when no ambiguity can occur. 
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Ma, Trudingcr and Wang's A3w condition and its strict version A3s then assert 

(1.3) A3w: cross(p, p) > for all p © p with h(p © p, p © p) = 0; 

(1.4) A3s: in addition, cross(p,p) = in p.3[) implies p = or p = 0. 

A3w and A3s are also called weak regularity and strict regularity, respectively. If 
(N, h) satisfies the inequality cross(p, p) > for all p(Bp (whether or not h(p(Bp, p(Bp) 
vanishes), then (N, h) is said to be non-negatively cross-curved. Loeper [Llj showed 
that A3w is necessary and A3s sufficient for continuity of optimal maps between 
suitable measures: without A3w there are discontinuous optimal maps between 
smooth measures p, p on nice domains. Trudinger and Wang [TWlj had already 
shown the sufficiency of A3w for continuity (indeed smooth differentiability) of op- 
timal mappings, under much stronger smoothness and convexity restrictions on p 
and p. These restrictions on p, p are relaxed in two dimensions by Figalli and Loeper 
|FL| . where they showed a continuity result of optimal maps under A3w. It still 
remains an open question to show such continuity result in higher dimensions. How- 
ever, in a separate work with Figalli [FKM , we show non-negative cross-curvature 
allows these smoothness and convexity restrictions to be relaxed without sacrificing 
Holder continuity of optimal maps and without assuming A3s, thus obtaining a 
continuity theory for smooth cost functions parallelling and extending Caffarelli's 
c(x,y) = \\x — y\ 2 result [Calj . 

For the Riemannian distance squared cost c(x,x) = dist 2 (x, x)/2 on a Rieman- 
nain manifold M = M, an isometric copy of M is embedded totally geodesically 
as the diagonal of M X M with respect to the pseudo-Riemannian metric h. Along 
this diagonal, 

4 

cross,^) = -Rm(p Ap,pAp) 

where Rm denotes the curvature operator of M (see jKMj for details) . This pro- 
vides some geometric intuition motivating Loeper 's result |L1| that A3w implies 
nonnegative sectional curvature of the Riemannian metric. (However, the nonncg- 
ative/positive curvature does not imply A3w, as shown by a counterexample in 
[K].) Loeper also verified A3s for the standard round sphere [L2] and used it to 
obtain C 1 ! max {M«-i} anc [ (j^a re g U i ar jty results for optimal maps in this spherical 
setting. 

1.2. Main results. Throughout this paper, if not specified, each Riemannian man- 
ifold M is assumed to be complete and to be equipped with the cost function 
c = ^dist 2 and this cost induces the pseudo metric h on N = M x M \ cut-locus. 
A Riemannian manifold M is said to be A3w / A3s / non-negatively cross-curved 
if (A, h) satisfies the corresponding cross-curvature condition. 

Our main results provide methods of generating new examples of Riemannian 
manifolds which satisfy non-negative cross-curvature and thus A3w. As announced 
in [KM) , we show: 

Theorem 1.2. (Products and submersions, particularly of round spheres) 

1) S n with its standard round metric is non-negatively cross-curved. 
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2) Let 7r : M — > M be a Riemannian submersion (see Definition ^- If M is A3w 
/ A3s / non-negatively cross-curved then so is M . 

3) For product M + x A/_ of Riemannian manifolds, if each factor M± is non- 
negatively cross-curved, then the resulting manifold M + x M- is non-negatively 
cross-curved, thus A3w holds (but never A3s). 

4-) Moreover, if either of the factors above fail to be non-negatively cross-curved 
then the product M+ x Af_ fails to be A3w. 

Proof. The proof of assertion 1) is a calculation given in Section [6l Assertion 2) 
is shown in Section [4] Assertion 3) and 4) are easy facts which are explained in 
Section [3] in detail. □ 

Remark 1.3. Following our announcement of Theorem 11.21 [KM . Figalli and Rif- 
ford |FR] gave an alternate proof of result 1) in a different and slightly stronger 
form. 

As a byproduct of our method, we obtain an O'Neill type inequality for cross- 
curvature in Riemannian submersions (see Theorem 14. 5|) . This verifies that Rie- 
mannian submersion quotients of the round sphere all satisfy A3s, from which the 
regularity of optimal maps follows on notable examples such as on complex projec- 
tive spaces CP™ with the Fubini-Study metric: see Section O This is a new result 
which is not covered by other discrete quotient cases of Delanoe and Ge |DG] or 
Figalli and Rifford [FR] , 

As another important consequence, the Riemannian product of round spheres 
and Euclidean space S" 1 x • • • x S' lfc x R' and its Riemannian submersion quo- 
tients, all satisfy cross-curvature non-negativity and thus A3w. Since this rules 
out the known counterexamples to regularity |L1| , optimism combines with a lack 
of imagination to lead to the conjecture that regularity of optimal mappings also 
holds in such settings. Together with the perturbations [DGj |LVj [FRj of the round 
sphere and its discrete quotients discussed by Delanoe and Ge, Loeper and Villani, 
and Figalli and Rifford, for which the continuity [FRj or regularity |LV] [DGj of 
optimal maps is already shown, these presently form the only examples of non-flat 
Riemannian manifolds which are known to be A3w. (In flat manifold case, A3w 
is trivial, and in fact the cross-curvature vanishes everywhere. Regularity of opti- 
mal maps in this flat case is known by Cordero-Erausquin [Co] applying Caffarelli's 
result [CaX][Ca2] .) 

As far as we know, it remains an open challenge to show regularity of opti- 
mal maps on the new tensor product type examples S™ 1 x • • ■ x S" fe x R' (when 
the supports of the source and target measure are the whole domain). However, 
Loeper's maximum principle (and a stronger convexity statement) is easily verified 
on these examples using the results and methods of [KMj : see Corollary 12.111 and 
Remark 13.41 We hope these key ingredients will make it possible to address the 
regularity issue in a subsequent work. 

This paper is organized as follows: Section [2] is devoted to preliminary notions 
and facts. Some important geometric implications of cross-curvature non-negativity 
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( Theorem 12.101 and Corollary 12. lip are shown also. In Section [3l the tensor prod- 
uct construction of costs is explained. Section [4] discusses the relation between 
Riemannian submersion and cross-curvature. Notably, we derive an O'Neil type 
inequality for cross-curvature (see Theorem 14. 5[) . Section [5] shows regularity of 
optimal maps on the complex projective spaces CP™. Section [6] shows the cross- 
curvature non-negativity of the standard round sphere. 

Acknowledgements. We thank Philippe Delanoe, Yuxin Ge, Alessio Figalli, Gregoire 
Loeper, Neil Trudinger, Cedric Villani and Xu-Jia Wang for useful discussions and 
timely exchanges of preprints. We are also grateful to the participants of Fields 
Analysis Working Group Seminar 2006-2008 for the stimulating research environ- 
ment they helped to create. 

2. Preliminaries 

In this section, we recall basic terminology and facts from |KM] (also |MTWj 
[TWlj |L1| ). Theorem 12. 101 and Corollary 12.111 are not stated there, but are easy 
consequences of the method in [KM] . to which we also refer the reader for proofs 
of the other results summarized below. For improvements to these results, see also 
[V2].Let JVcMxMbean open set where the cost function c G C 4 (iV). 

Definition 2.1 (visible sets, twist condition and non-degeneracy). We de- 
fine the visible sets: 

N(x) = {xe M | (x,x) G N}, 
N(x) = {x G M | (x,x) G N}. 

Throughout this section it is assumed as in [MTWj that for all (x,x) G N, 

Al Dc(x, ■) : N(x) -> T*M, Dc(-,x) : N(x) -> T*M are injective; 
A2 DDc is non-degenerate. 

We also call Al the (bi-)twist condition and A2 non- degeneracy. 

We recall an important map of Ma, Trudinger & Wang [MTWj , called the cost- 
exponential by Loeper |L1| , which coincides with the Riemannian exponential map 
for c = ^dist 2 . 

Definition 2.2. (cost exponential) If c G C 2 (N) is twisted (Al), we define the 
c- exponential on 

Don^c-ExpJ := - Dc(x, N(x)) 
(2.1) ={p* e T*M | p* = -Dc(x, x) for some x G jV(ir)} 

by c-Exp x p* = x if p* = —Dc(x,x). Non-degeneracy (A2) then implies the c- 
exponential is a diffeomorphism from Dom(c-Exp x ) C T*M onto N(x) C M. 
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Remark 2.3 (differential of c-Exp). Linearizing the cost exponential x = 
c-Exp^p* around p* G T*M we obtain a map c-Exp^^,, : T*M — > T S M given 
explicitly by 

(2.2) oExp (S)S(t)) ,(p*(t)) = x{t) forp*(t) = -Dcfosft)). 

Equivalently c-Exp^ s ^ = —DDc(x, a;) -1 , meaning the inverse tensor ^ft _1 to the 
metric — which gives the pseudo-Riemannian correspondence between the 

tangent and cotangent spaces to N C M X M — also carries covectors forward 
through the cost-exponential. 

The following lemma characterizes Ma, Trudinger & Wang's c-segments as geodesies 
of ft. 

Lemma 2.4. (the c-segments of [MTWj are geodesies) Use a twisted (Al ) and 
non-degenerate (A.2) cost c G C 4 (iV) to define a pseudo-metric on the domain 
N C M X M . Fix x G M . For each line segment (1 — s)p* + sq* € Dom(c-Exp x ), 
s G [0,1], the curve 

s e [0, 1] — ► a(s) := (x,c-Exp x ((l - s)p* + sq*)) 

is an affinely parameterized null geodesic in (N,h). Conversely, every geodesic seg- 
ment in the totally geodesic submanifold {x} x N(x) can be parameterized locally 
in this way. 

To see some relevant geodesic equations in local coordinates, given s G [0,1], 
introduce coordinates on M and M around o~(sq) so that nearby, the curve <r(s) 
can be represented in the form (x 1 , . . . , x n ,x 1 (s), . . . , x n (sj). Differentiating the 
definition of the cost exponential 

(2.3) 0=(l-s)p*+sq*+ Ci (a(s)) 
twice with respect to s yields 

(2.4) = djx k + c m x j x k 

for each i = 1, . . . , n. This equation will be used later. 

Regarding curvature of the metric ft, the following fact is fundamental. 

Lemma 2.5. (Non-tensorial expression for curvature) Use a non- degenerate 
cost c £ C A {N) to define a pseudo-metric (jl.ip on the domain N C M X M. Let 
(s,t) G [— 1,1] 2 — ► (x(s), x(t)) € N be a surface containing two curves o~(s) = 
(x(s),x(0)) and r(t) = (x(0),x(t)) through (x(0),x(0)). Then © x(0) defines a 
parallel vector-field along o~(s). Moreover, if s G [—1, 1] — ► o~(s) G N is a geodesic 
in (iV, ft) then 

d 4 

c(x(s),x(t)) = cross( x ( ), S (o))(i;(0),x(0)). 



(2.5) - 2 



ds 2 dt 2 



Note in this lemma that only one curve a(s) needs to be geodesic in (TV, ft). As a 
consequence of this result, the conditions A3w/s can alternately be characterized 



CURVATURE OF OPTIMAL TRANSPORT 



7 



by the concavity / 2-uniform concavity for each x G M and q$ G Dom(c-Exp a; ) of 
the function 

(2.6) q* G T* X M — pV' C!J (x, c-Exp x (g* + g*)). 

restricted to g* in the nullspace of p G T X M . Cross-curvature nonnegativity asserts 
this concavity extends to all q* (not necessarily in the nullspace of p) such that 
<7o + Q* G Dom(c-Exp x ). 

Before recalling the important geometric implications of the curvature properties 
of h, let us define: 

Definition 2.6. (Illuminated set) Given (x,x) G N, let V(x,x) C M denote 
those points y G N(x) for which there exists a geodesic curve from (x, x) to (y, x) 
in N(x) x {x}. 

We now state a version of Loeper's maximum principle. Although the state- 
ments below are from |KMj , where a direct elementary proof is given using pseudo- 
Riemannian geometry, the fundamental form of the equivalence (Theorem [2TTJ) be- 
low was deduced by Loeper |Llj in the simpler setting N = M x M C R n x R" from 
results of Trudinger and Wang TWT] |TW2] . We visualize his maximum principle 
as asserting that the double-mountain max[/o, /i] stays above the sliding mountain 
ft{y) '■= —c(y, x(t)) + c(x, x(t)), hence refer to it by the acronym DASM. 

Theorem 2.7 (A3w local DASM). Let h be the pseudo-Riemannian metric 
on N C M X M induced from the non-degenerate cost c G C 4 (A) as in (jl.ip . The 
following are equivalent. 

1. (N,h) satisfies A3w. 

2. (local DASM) For any h-geodesic a : t G [0, 1] — ► (x,x(t)) G and suffi- 
ciently small neighbourhood U C V a of x, where V a := Ho<t<iV(x, x(t)) is from 
Definition II?. 61 the sliding mountain ft(y) '■— —c(y,x(t)) + c(x,x(t)) satisfies the 
maximum principle 

ft(y) <max[/ ,/i](y) for0<t<l. 
Then one can prove as in K M the following. 

Theorem 2.8 (A3w + c-convexity of domains => DASM). Let h be the 

pseudo-Riemannian metric on N C M x M induced by the non- degenerate cost 
c G C 4 (N) n C(M x M) as in JlTJ. Suppose (N,h) is A3w, and the set V a := 
^0<t<iV(x,x(t)) from Definition \2.6] is dense in M for some h-geodesic a : t G 
[0, 1] — ► (x, x(t)) G N. For any y in M , the sliding mountain ft(y) '■= —c(y, x(t)) + 
c(x,x(t)) satisfies the maximum principle 

(2.7) DASM: f t (y) < max[/ , fi](y) for < t < 1. 

Remark 2.9 (Relating density of the illuminated set to c-convexity of the 
domain). In case N = M x M, the density hypothesis of the preceding theorem 
holds whenever Dc(M,x) C T^M is convex for each x G M, since then V a = M. 
This is the case considered initially by Trudinger and Wang [TWlj and Loeper [Ll] . 
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For our argument, it is enough that Dc{M, x(t)) be star-shaped around Dc(x, x(t)) 
for each t £ [0, 1]. Conversely, V(x, x) = M for each x £ M implies convexity of 
Dc(M, x). 

Our pseudo-Riemannian method |KM| makes it equally possible to deduce fur- 
ther geometric implications of the cross-curvature condition, as the next theorem 
and corollary show. 

Theorem 2.10 (nonnegative cross-curvature local time-convex sliding 
mountain). Let h be the pseudo-Riemannian metric on N C M X M induced from 
the non-degenerate cost c G C A {N) as in The following are equivalent. 

1. (N,h) is non-negatively cross-curved. 

2. For each h-geodesic a : t £ [0, 1] ► (x,x(t)) £ N and sufficiently small neigh- 
borhood U C V a of x £ U, where V a :— r\o<t<lV(x,x(t)) is from Definition \2.b\ the 
sliding mountain ft{y) '■= — c(y, x(t)) + c(x, x(t)) is a convex function of t £ [0,1] 
for each y£U, i.e., J^/ t (j/) > for < t < 1. 

Proof: (1 implies 2). Fix an arbitrary /i-geodesic t : [0,1] — > (x,x(t)) £ N. Let 
U be chosen open so that x £ U C V a . The existence of such a U is elementary, 
tedious, and independent of hypothesis 1: it requires checking — at least for 5 
sufficiently small depending on T £ [0, 1] — that dT-5<t<T+8V{x, x{t)) contains a 
neighbourhood of x, as we now do. Choose coordinates on M near x. Since the 
geodesic a is compact in the open set N, some coordinate ball satisfies B r (x) x 
{x(t)} C N for all t £ [0, 1]. Since the coordinate charts x £ M — ► Dc(x,x(t)) £ 
T* {t) M are C 2 smooth functions of (x, t) £ M x [0, 1], taking r = r(T) and 5(T) > 
sufficiently small ensure Dc(B r (x), x(t)) is convex for each t within 6(T) of T. 
This convexity implies B r {x) C V(x,x(t)). Extracting a finite subcover ufl 1 (Tj — 
5(Ti), Tj + 6(Tij) of [0, 1] and taking U = r\i< N B T{Ti) {x) will suffice. 

Now define /<(•) := — c(-,x(t)) + c(x,x(t)). Fix arbitrary to £ [0, 1] and y £ U. 
The fact U C V a guarantees an /i-geodesic s : [0,1] — > (x(s),x(to)) € with 
x(0) = x and x(l) = y. Define an auxiliary function g(s) :— ^ \ t _ to ft(x(s)), 
which shall be shown to be non-negative for s £ [0,1]. By Lemma [2.51 and the 
cross-curvature non-negativity, 

d 2 a 

(2.8) -4f >0, for s e [0,1]. 

In particular, g(s) is convex. It is clear that g(0) — ^r| t _ to /t(x) = 0. We also 
claim g'(0) = 0: introducing coordinates x l ,...,x n around x — x(0) on M and 
x , . . . , x n around x(to) on M, we compute 

= -{c- ik (x(0),x{t o ))x l + c l]k {x(0),x{t ))x l x : >)x k 

s=0 

= 



(2.9) ^ 
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by the h-geodesic equation ([2^4]) for t £ [0,1] -> (x,x(t)) £ N. From (|2~8| . this 
shows g(s) > 0, for Q < s < 1, in particular at s = 1, 



2 



3 

9W = aj? My) 



> o. 

t=t 



Since we have used U d V a but not the openness or smallness of t/, we have actually 
deduced convexity of t 6 [0, 1] — ► /t(j/) for all y £ V CT . A fortiori, 1 => 2. A 

Proof: (2 implies 1). Fix (x, x) £ N, p © p £ Tr x S .)N, It shall be shown that 
sec(sc,s)({p 0) A (Offip)) > 0. Choose an /i-geodesic t £ [-1,1] — * (a;,x(t)) € iV, 
with x(0) = x, x(0) = p, and a curve s S [—1, 1] — » (y(s),x) £ N, with y(0) = x, 
y(0) = p. Let /(•) := — c(-, + c(x, Suppose 

(2-10) g(s) := ^ My(s)) > 0, 

for s £ [—1,1]: this holds if the property 2 is assumed and the curve y(s) is 

chosen inside the neighborhood U of x constructed at the outset. Note g(0) = 
$i\ t=0 ft{x) = 0. Thus, from (LTTOll . g'(0) = and 



d 



9(s) 

s=0 

cross (Xj2) (y(0),i(0)). 



2 

The last equality comes from Lemma 12.51 This completes the proof 2 1 . □ 

Corollary 2.11 (non-negative cross- curvature + c-convexity of domains 
time-convex sliding mountain). Suppose the pseudo-Riemannian metric h 
induced by the non- degenerate cost c £ C 4 {N)nC{M x M) on N C M x M in (fTTT]) 
is non-negatively cross-curved. If for some h-geodesic a : t € [0, 1] — ► (x, i(<)) £ N 
the set V a := rio<t<i V(x, x(t)) /rom Definition \2.b\ is dense in M ; then for each 
y £ M, the sliding mountain ft{y) '■= —c(y,x(t)) + c(x,x(t)) satisfies 

(2.11) ft(y)<(l-t)fo(y)+tf 1 (y) for0<f<l. 

Proof. For y £ V a , inequality (|2. 1 1[) was established while deriving implication 1 
=>■ 2 of Theorem 12. 101 The inequality extends to y £ M by the density of V a and 
the continuity of c on M x M. □ 

Remark 2.12. The density condition of Theorem 12.81 and Corollary 12.111 can be 

further weakened by the works of Figalli, Loeper and Villani |FV| |LV] (see |V2[ 
Theorem 12.36]). 



3. Tensor products of pseudo-Riemannian metrics 

This section contains the proofs of assertions 3) and 4) in Theorem ll.2l Assertion 
3) is actually a corollary of a more general theorem, which does not require the 
transportation cost c or manifold M to be Riemannian: 



10 



YOUNG-HEON KIM AND ROBERT J. MCCANN 



Theorem 3.1. (Products preserve non-negative cross-curvature) Let c± £ 

C (N±) be non-degenerate non- negatively cross-curved costs on two manifolds N± C 
M± x M±. Then c(x + , x_, x + , x_) = c + (x + ,x + ) + c_(x_,x_) is non- degenerate 
and non-negatively cross-curved on (x+, x+, x_, x_) £ N + x 7V_, but never A3s. 

Proof. If the cost functions c± : N± C M± x M± — > R define non-degenerate 
metrics h±, then so does the cost c(x+, x_, x+, xJ) = c+(x+,x+) + c_(x_,x_) 
on N = N + x N- , since the corresponding metric h separates into block diagonal 
components h± with non- vanishing determinants. In this sense the geometry (iV, h) 
is the pseudo-Riemannian tensor product of the geometries (N±, h±). It follows that 
the product (7+(s),7_(s)) of geodesies s £ [—1,1] — > 7±( s ) m N± is a geodesic in 
N = N + x AT_ . Lemma 12.51 then implies 



(3.1) cross(p+ ®p-,p + = cioss + (p + ,p + ) + cross_ (p_ , p_ ) 

for all (x + ,x+,x_,x_) € N and ( P+ ffip_)ffl(p + fflp_) e r (a , +i 3 +iai _ )S _)JV". The 
proof is completed by observing that non-negativity of both summands guarantees 
the same for their sum. Choosing p_ = p + = in identity (|3.1[) demonstrates that 
the tensor product cost c cannot satisfy A3s. □ 

Proof of Theorem 11.21( 3). Given two Riemannian manifolds M± = M±, de- 
note their geodesic distances squared by c±(x±,x±) — hdiat±(x±, x±) respectively, 
and let c(x+, x~, x+, x_) — ^dist 2 ((x+, x_), (x+, x_)) denote the geodesic dis- 
tance squared on M + x M_ equipped with the Riemannian product metric. Then 
c(x+, X—, x+, x_) = c + (x + , x + ) + c_(x_ , x_), and the non-negative cross-curvature 
of c on ]V = Ar + x iV_ follows from Theorem 13. 1( taking N± to be the domains 
where c± are smooth. This concludes the proof of Theorem [L~2^ 3). □ 

The assertion 4 in Theorem 11.21 follows easily from the next lemma. 

Lemma 3.2 (Certain cross-curvatures vanish in any Riemannian setting). 

Define the pseudo-metric (|1.1[) using the cost c = rf 2 /2 on a Riemannian manifold 
(M, d). Each point (x, x)eN = MxM\ cut-locus with x ^ x admits a 2-plane 
whose cross- curvature vanishes. For example, letting t £ [0,1] — > 7(4) £ -A/ be a 
geodesic from x — 7(0) to x — 7(1) yields cross(7(0), 7(1)) = but 

fe(7(0)ffi7(l),7(0)ffi7(l))-^,S) 2 . 



Proof. Any (afiinely parameterized) Riemannian geodesic satisfies d(-j(s), "f(t)) 
\s-t\d(j 
compute 



|s-t|d(7(0),7(l)). Defining f(s,t) = c{^{s),j(t)) = \s - t\ 2 d(x, x) 2 /2 allows us to 



r\2 £ 

(3.2) ft( 7 (0)®7(l),7(0)®7(l)) 



dsdt 



(s,t)=(0,l) 



immediately. Moreover, the fact that t £ [0, 1] — > (x,j(t)) is an ^.-geodesic yields 
the vanishing of cross(7(0), 7(1)) = — 29 4 //9s 2 9i 2 |( s t ) = (o.i) via Lemma |2~51 □ 
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Proof of Theorem [TT2l( 4). Without loss of generality, assume M+ fails to be 
non-negatively cross-curved. Then cross+ (p+ , p+ ) < for some (p+ , p+ ) G T( x+ <m \N+. 
Noting 

h +(p+ © ±P+,P+ ©±P+) = ±h + {p + ®p+,p+ ®p+) 

but cross+ (p+ , ±p + ) = cross+ (p+ , p+ ) , we may assume that h+ (p+ ® p+ , p+ ©P+ ) < 
0. Pick any nontrivial geodesic 7_ in M_. Noting /i_(7_(0) © 7 _(1), 7 _(0) © 
7-(l)) = rf 2 (7(0),7(l)) > from ([372]) . one can choose A G R so that 

/i((p + ffi A 7 _(0)) © (p + © A 7 _(l)), (p+ © A 7 _(0)) © (p + © A 7 _(l))) 
= h+(p+ ®p+,p+ ®p+) + A 2 /i_( 7 _(0) © 7_(1),7_(0) © 7 _(1)) 
= 0. 



However, from (|3.1[) and Lemma 13 . 21 

cross(p + © A7_(0),p+ © A 7 _(l)) 
= cross+(p + ,p+) + A 4 cross_(7_(0),7_(l)) 
= cross + (p + ,p + ) < 0. 

This completes the proof that (N + x JV_, h + @h-) fails to be weakly regular A3w 
unless both (N±, h±) are non-negatively cross-curved. □ 

At present, the authors do not know any Riemannian manifold which is A3w yet 
fails to be non-negatively cross-curved. On the other hand, the proof of Th.eorem ll.2l 
(4) adapts to non-Riemannian cost functions as in the next example, where it shows 
that the tensor product of two (or more) costs on manifolds which are not all non- 
negatively cross-curved is surely not A3w. 

Example 3.3 (A3s x A3s =f> A3w). Let c±(x,x) = —log\x—$\ onM±xM±\A, 
M± = R", A := {(x,x) | (x, x) G M± x M±}. This logarithmic cost function is 
known to be A3s |MTW| [TWlj but it does not induce non-negatively cross-curved 
pseudo-metric as indicated, e.g. in [KM, Example 3.5]. To see this fact one uses 
(|2.6[) . whose righthand-side coincides for this logarithmic cost with 



, n „ u , =2((q + q)*P 1 ) 2 - \p\ 2 \{q + q)*\ 2 



where f(x — x) := —log\x — x\. Here the righthand-side is strictly convex with 
respect to q* along the Euclidean line parallel to p, but strictly concave along the 
nullspace of p. As a result cross±(p±, Xp±) < for p± parallel to DDc(x, x)p± as 
Euclidean vectors and A ^ 0. On the other hand, 

h+{p+ ®P+,P+ ®P+) + h^(p^ © Xp-,p- © Xp-) 



vanishes for some non-zero A G R. From (|3.ip . the pseudo-Ricmaniann metric h 
induced by the cost c((x, x), (y, y)) = c + (x + ,x + ) + c_(x_,a;_) on (M + x M_) x 
(M + x M_) then fails to be A3w. 
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Remark 3.4 (tensor product examples, Loeper's maximum principle, and 
time-convexity of the sliding mountain). As mentioned in the introduction, 
Loeper's maximum principle (DASM) is a key property for the regularity of op- 
timal maps. The conclusions of Theorem EH and Corollary I2TT1 (hence DASM) 
hold for the distance squared cost on the Riemannian product of round spheres 
M = M = S" 1 x • • • x S" fc x R 1 — thus also on its Riemannian submersions (see 
Theorem 14. 8|) . To see this, first note that by the result of present section and 
Section [51 M satisfies non-negative cross-curvature on N = M x M \ cut-locus. 
The density condition of rio<t<i V(x, x(t)) is easily checked since the cut locus of 
one point in this example is a smooth submanifold of codimension greater than 
or equal to 2. This new global result illustrates an advantage of our method over 
other approaches [TWT] [TW2] [Ll] . where one would require a regularity result for 
optimal maps (or some a priori estimates) to obtain the conclusion of Theorem EU 
For example, to implement these other approaches for the manifolds of this ten- 
sor product example, one would need to establish that an optimal map remains 
uniformly away from the cut locus, as is currently known only for a single sphere 
M = M — S™ from work of Delanoe & Loeper [PL] (alternately [Llj or Appendix 
of [KM]), for the case of perturbations |DG] [LV] |V3j . and for some Riemannian 
submersion quotients (see Section[5]) of S", including the complex projective spaces 
CP™ with Fubini-Study metric. To the best of our knowledge, no one has yet suc- 
ceeded in establishing regularity results for this tensor product example, though as 
mentioned above we hope to resolve this in a subsequent work. 

4. Riemannian submersions and cross-curvature 

In this section we prove the assertion 2) in Theorem 11.21 The key result is an 
O'Neill type inequality contained in Theorem 14. 51 which compares cross-curvatures 
in Riemannian submersion. Thcorem l4.8l deals with the survival of global properties 
such as Loeper's maximum principle (DASM) and time-convexity of the sliding 
mountain, under Riemannian submersion. These results are applied to show the 
regularity of optimal maps on the Riemannian submersion quotients of the round 
sphere, for example, the complex projective space CP™ with Fubini-Study metric 
(see Section [5]). From now on we focus on the case of Riemannian manifolds (with 
c = idist 2 ). In this case the c-exponential map coincides with the Riemannian 
exponential map: c-Exp = exp. 

Recall the definition and basic facts of Riemannian submersion. 

Definition 4.1 (Riemannian submersion). (See |CEJ.) A surjective differen- 
tiable map 7r : M — > B from a Riemannian manifold M onto a Riemannian manifold 
B is said to be a Riemannian submersion if the following hold: 

• 7r is a submersion, i.e. dn : T X M — » T w r x \B is surjective for each x £ M ; 

• for the orthogonal decomposition T X M = ker dirQ) (ker dir) 1 - for each x £ M, 
dn\ ( kcrd7t .)j_ is an isometry. 

The subspaces V := kere?7r, Ti := {kendir) 1 ' are called vertical and horizontal sub- 
spaces, respectively. For each v £ T^B, b £ B, there exists its unique horizontal 
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lift v E T X M n TL for each x E 7r _1 (6) such that dnr(v) — v. We use the metric 
identifications T b *B — T b B and T*M = T X M to extend the definition of a horizon- 
tal lift to cotangent vectors. For each piecewise smooth curve 7 : [0, 1] — > B and 
x € 7r _1 (7(0)), there exists its horizontal lift 7 : [0, 1] — ► M such that 7(0) = x, 
tt(j) = 7, 7(i) G 7i. Moreover, 7 is a geodesic if and only if its horizontal lift 
7 is. If 7 is minimal, then so is 7. This property yields, for the horizontal lifts 
v E T X M (7 H and x E tt" 1 (b) of v E T b B and b E B, 

(4.1) ir(exp x v) = exp b v. 
Regarding Riemannian distance, 

(4.2) <3M M (x,y) > dist s (7r(ar), n(y)). 

We call M the total space of the Riemannian submersion and B the base of the 
Riemannian submersion or the Riemannian submersion quotient. 

Many examples of Riemannian submersions may be found in Cheeger & Ebin [CEj 
or Bess [Be]. Every Riemannian covering projection is obviously a Riemannian 
submersion. Other important examples are Hopf fibrations such as complex and 
quaternionic projective spaces CP" 1 and HP™, where the standard round sphere 
S n (sectional curvature = 1) is the total space. 

Example 4.2 (Hopf fibrations). (See [Be] pages 257-258.) it : S 2m+1 -> CP™, 
7T : 5' 4m + 3 -> HP m . The base spaces CP'" and HP m have real dimensions 2m, 
4m, respectively, and have non-isotropic sectional curvatures K, 1 < K < 4. 

Our goal in this section is to compare the cross-curvature of the base space with 
that of the total space of Riemannian submersion. For this purpose, we use the 
definition above to assign to each pair of points in the base space a corresponding 
pair of points horizontally lifted to the total space. 

Definition 4.3 (horizontal lift of a pair of points). Let ir : M — > B be a 

Riemannian submersion. For each pair of points (x, x) E Nb = B x B \ cut-locus, 
the unique geodesic 7 : [0, 1] — > B with 7(0) = x, 7(1) = x has its horizontal lift 7 
with 7(0) = x, 7(1) = x. Then the pair (x, x) E Nm = M x M \ cut-locus is said 
to be a horizontal lift of (x,x). 

To have a bit more general conclusion in the next theorem, we observe as in [M] 

Lemma 4.4 (radial cost-exponential), [see [Ml Theorem 13]/ Let f : R — > R 
be a strictly convex smooth function, thus the derivative f has an inverse function. 
Further assume that f is strictly increasing on R + . Let M — M be a Riemannian 
manifold, and c{x,x) = /(dist(a;, x)) a cost function defined on M x M, where dist 
denotes the Riemannian distance. The c-Exp and exp have the relation 

fUTH\p\). 



(4.3) c-Ex Px {p* 




-p) for 7^ p* E Dom c-Exp x , 
for p* = 

where the vector p and the co-vector p* are identified by the Riemannian metric. 
This cost c induces a pseudo-Riemannian metric hu on Nm '■= M X M \ cut-locus 
as in (jl.ip which is A2 non-degenerate and Al bi-twisted. 
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Proof. For f{p) — p , this fact is well-known. For general /, it follows from 

/ / (dist(x,x)) n , dist 2 _ 
- DC{X ' X} = ~ dist(x,x) D ( — ^ x ^ 

where the inverse exists for strictly convex function /. The last assertion 

can be checked directly from Definition 12.11 This completes the proof. □ 



Now we are ready to show our main theorem in this section. 

Theorem 4.5 (cross-curvature and Riemannian submersion). Let n : M — > 

B be a Riemannian submersion from M to B. Let f : R — > R be a strictly 
convex smooth function that is strictly increasing on R + , and cm = f o distM and 
cb = f ° dists be cost functions defined on M x M and B x B, where dist denotes 
the Riemannian distance of the corresponding manifold; cm cind cb induce the 
pseudo- Riemannian metrics Km and hs on Nm :— M x M \ cut-locus and Nb ■= 
B x B\ cut-locus respectively, as in Fix (x,x) £ Nb and let (x,x) £ Nm be 

a horizontal lift of (x,x). Given v v € T^ x S )Nb there exists w w £ T^^Nm 
with diTx(w) = v, dir^(w) = v, such that 

(4.4) h B (v®v,v®v) = h M (w®w,w® w) 
and 

(4.5) crossg^Ks) > crossg^^). 

For example, it suffices to take w* = —DDcm(x, x)w £ T-M to be the horizontal 
lift of v* = —DDcb{x,x)v £ T*B and w* = —DDcm(x,x)w £ TSM to be the 
horizontal lift of v* = -DDcb(x,x)v £ T?B. 

Proof of (|4.5[) . Let (x, x) be a horizontal lift of (x, x) £ Nb, and define q* and q* by 
cs-Exp^g* — x and cs-Exp s g* = x. Then dists(a;,a;) = distj\/ (£, x) and it follows 
from (|4.3| that the horizontal lifts q* of q* and q* of q* satisfy CAf-Exp^g* = x 
and CM-Exp=q* = x. To fixed v £ T X B and v £ T S B correspond v* £ T^B and 
v* £ T*B such that v* = —DDcb(x,x)v and v* = —DDcb{x,x)v. Equivalently, 
v* © v* = 2h B (v <Sv, ■). 

Now let £ : (s,t) £ [—1, l] 2 — > (x(s),x(t)) £ Nb be the surface given by x(s) — 
cg-Exp 2 {q* + sv*) andx(t) — cs-Exp 2 .(g* +tv*). By Lemma [241 the curves a(s) — 
(x(s),x) and r(t) — (x,x(t)) through (x(0), x(0j) = (x,x) are /i^-geodesics. Lift £ 
to £ : (s, t) £ [—1, l] 2 — > {i{s), x(t)) £ Nm by setting x(s) = CA/-Exp=(g* + sw*), 
and x(t) — CM-Expj(g* + tw*), where w* and w* are the horizontal lifts of v* 
and v* respectively. Thus, &(s) — (x(s),x) and f(t) — (x,x(t)) are /iM-geodesics, 
with cr(0) = (x,x) = t(0). Moreover, n(x(s)) = x(s) and ir(x(t)) = x(t) from 
(|4.3[) and (14. ip . so taking w :— i(0) and w := x(Q) yields d-Kx(w) = i(0) = v and 
dir^(w) = x(0) = v. Notice that 

-DDc B (x(0),x)x(0) =«*, -S£)c B (x,x(0))x(0) = u* ; 

-DSc M (x(0),S)^(0) = w*, -L>Z)cAf(x,l(0))i(0) = w*. 
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Define an auxiliary function 

F(a,t) := c M (x(s),x(t)) - c B (x(s),x(t)). 

From Lemma [2?5| the inequality d ®i gt -i F(s, t)\, Q QS) > will imply (|4.5|) . First, 
observe from (|4.2p and the monotonicity of / that cb(x(s), x(t)) < CM(x(s),x(t)), 
thus, 

(4.6) F(s,t)>0 for (s, t) e [-1, l] 2 . 

We shall verify the desired inequalities by computing the Taylor expansion of F at 
(0,0) to fourth order. First observe that F(s,0) = = F(0,i) for all |s|, |t| < 1, so 

F(s, t) = fast + f 2 is 2 t + / 12 si 2 + f 31 s 3 t + f 22 s 2 t 2 + f 13 st 3 + 0(\s\ + \t\) 5 

as \s\ + \t\ — > 0. Since F(s, ±s) > for |s| < 1, we deduce the vanishing of /n, /12 
and /21 in turn, and the inequalities /31 + f 22 + /13 > and ^22 ~ /31 — /13 > 0. 
This implies /22 > as desired. (Although not needed here, /31 = follows from 
F(s, s 2 ) > 0, and /13 vanishes similarly.) Noting 

d 2 F 

= i(0)DDc M (x, x)x(0) - x(0)DDc B {x, x)x(0) 
= -h M (x(Q) © 1(0), i(0) © 1(0)) + h B (x(0) © x(0),x(0) © x(0)) 
we have established f|4.4[) en passant to complete the proof. □ 



Before stating a corollary of this theorem, let's make a provisional definition 
which can serve as a strict cross-curvature condition for a Riemannian manifold. 
Notice from Lemma l3.2l that for each pair of points in a Riemannian manifold there 
are tangent vectors with zero cross-curvature. 

Definition 4.6 (almost positive cross-curvature). A Riemannian manifold M 
with positive sectional curvature is said to be almost positively cross-curved if for 
each (x, x)EN = MxM \ cut-locus such that x ^ x and p(Bp £ T/ Xt g)N, 

(4.7) cross(p,p) > 

and the equality holds if and only if p and p are parallel to the velocity vectors 7(0) 
and 7(1), respectively, for the unique geodesic t E [0, 1] — > 7(f) G M from x to x. 

For example, the standard round sphere is almost positively cross-curved as shown 
in Section [5] 

Corollary 4.7 (A3w/A3s, non- negative/almost positive cross-curvature 
survive Riemannian submersion). Let it : M — * B be a Riemannian submer- 
sion. If the distance squared cost on M satisfies A3w, A3s, non-negative cross- 
curvature, or almost positive cross- curvature condition, then B satisfies the same 
condition, respectively. 
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Proof. The relevant inequalities for the cross-curvature follow directly from (14. 5|) . 
Let's consider especially the equality case of almost positive cross-curvature. As- 
sume M is almost positively cross-curved. Suppose 

crossg^ 3 = 

for x 7^ x, (x,x) S Nb, p®pG T/ x ^\Nb- Lift the unique minimizing geodesic 7 
linking 7(0) = x to 7(1) = itoa horizontal geodesic 7 joining x — 7(0) to x = 7(1). 
There is a unique choice p (Bp € T^^Nm such that each component p* and p* 
of p* (Bp* = h,M(p®P, •) is the horizontal lift of the corresponding component of 
p* (Bp* = fiB (p (B p, ■ ) ■ The preceding theorem asserts dir % (p) = p and dir^. (p) = p. 
Apply (|4.5p to get cross|^ r |^' ?iM '(p,p) = 0. Then by almost positive cross-curvedness 

of M, the vectors p and p are parallel to the velocities 7(0) and 7(1), respectively. 
This implies that the vector projections p and p are parallel to 7(0) = ^5(7(0)) 
and 7(1) = d7T|(7(l)), respectively, to complete the proof. □ 

Let's now turn to more global aspects of the distance squared cost function un- 
der Riemannian submersion. Though local DASM/local time-convex sliding 
mountain are equivalent to A3w/nonnegative cross-curvature, respectively (see 
Theorem 12 . 71 and Theorem l2.10|) , their global counterparts DASM/time-convex 
sliding mountain require additional conditions on the geometry of the domain 
(see Theorem 12.81 and Corollary 12. lip . The following theorem, uses a simple com- 
parison of distance to give a direct proof that both Loeper's maximum principle 
and time-convexity of the sliding mountain survive Riemannian submersion even 
in the absence of restrictions on domain geometry. 

Theorem 4.8 (Loeper's maximum principle and time-convexity of the 
sliding mountain survive Riemannian submersion). Let ir : M — > B be a 

Riemannian submersion. Compose a smooth strictly convex function f : R > R 

that is strictly increasing on R +; with the Riemannian distance on M to define 
a cost function cm = f distM- It induces a pseudo- Riemannian metric Km on 
Nm '■— M X M \ cut-locus as in (|1.1|) which is A2 non- degenerate and Al bi- 
twisted. Similarly cb = f ° diste defines a non-degenerate and bi-twisted pseudo- 
metric Kb on Nb := B x B \ cut-locus. Suppose that for each h,M-geodesic of the 
form t e [0, 1] — > (x,x(t)) in N M , ft(y) = -cm(v, x{t)) + c(x, x(t)) satisfies (|2.7|) 
(or (|2.11[) respectively) for each y € M. Then for each IiB-geodesic t € [0, 1] — ► 
(x,x(t)) € Nb, ft(y) = ~ c B(y,x(t)) + cb{x, x{t)) satisfies the same inequality for 
each y € B. 

Proof. Let t G [0, 1] — ► a(t) = (x, x(t j) 6 Nb be an /ig-geodesic and set x := x(0). 
Define the sliding mountain f t (-) :— —Cb(-, x(t))+CB(x, x(t)) on B. Identify tangent 
vectors with co-tangent vectors by the Riemannian metric. By Lemma 12.41 there 
exist p,q € T X B such that x(t) = c-Exp^p + tq). Lift p, q to horizontal vectors 
p, q at x S tt^ 1 {x). Let x{t) — c-Exp^^ + tq). From the Riemannian submersion 
property and Lemma 14.41 

c B (x,x(t)) = {f)- 1 (\p + tq\ B ) = c M (x,x(t)) for all te [0,1]; 
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and (x,x(t)) G Nm for t G [0,1]. This last property comes from the fact that 
Riemannian submersions lift the minimal geodesic from x to x(t) to the minimal 
geodesic from x to x(t). Thus, t G [0, 1] — > (x,x(tj) G AT M gives an /iM-geodesic. 
Define / t (-) = -c(-,z(*)) + c(5, §(*)). Fix i G [0, 1], y £ B. Let 7 : [0, 1] -> £ be a 
geodesic from 7(0) = x(t) to 7(1) = y. Let 7 be the horizontal lift of 7 such that 
7(0) = x(t). Let y := 7(1) G 7r _1 (y). Notice that 

c B (y,x(t)) = c M {y,x(t)); c B (y,x(s)) < c M (y,x(s)) for all s G [0,1]. 

The last inequality is from (|4.2[) and the monotonicity of /. Therefore, 

(4.8) / t (y) = f t (y); f s (y)<f s (y) for all a e [0, 1] . 

Now, assume s G [0, 1] — ► / s (y) is convex. Choosing s = 0, 1, from (|4.8[) . 

ft(y) = ft(y)<(i-t)f (y)+tf 1 (y) 
< (i- 0/o(i/) 

Since f G [0, 1] was arbitrary, the same convexity holds for t G [0, 1] — > /*(?/)• The 
survival of Loeper's maximum principle (DASM) follows by a similar argument. 

□ 



5. Regularity of optimal maps on the complex projective space CP" 

In this section, we establish regularity of optimal maps on the complex projective 
space CP™, in fact more generally on the Riemannian submersion quotient of the 
round sphere with strictly convex domain of exponential map. 

From Corollary 14.71 and Theorem 14.81 any Riemannian submersion quotient B 
of the round sphere satisfies, in particular, A3s and Loeper's maximum principle 
(DASM). Let's further assume that the domain of exponential map Dom(Exp :c ) 
is strictly convex for each x G B. Here the domain of exponential map is the 
special case of (|2.ip with c = idist 2 . We also identify vectors and co- vectors by 
the Riemannian metric. The Riemannian manifolds CP™ and HP™ in Example l4.2l 
all satisfy this condition. In these cases, the domain of the exponential map is the 
ball of radius ^ in the tangent space, and the conjugate locus coincides with the 
cut locus. Thus these manifolds are focal, in contrast to the nonfocal manifolds 
analyzed by Loeper and Villani in the preprint version of [LVj Pl Exactly the same 
method as in Appendix E of [KMj . which uses Loeper's argument jLlj . shows the 
Holder continuity of optimal maps for the source and target measures p, p on B, 
assuming there exists A > such that 

A<-^- and ^<1/A, 
a vol a vol 



^Remark added in revision: Shortly after we communicated the present manuscript to Loeper 
and Villani, we learned they had revised |LV) to address A3s manifolds which are not purely- 
focal. Our results of the preceding section establish CP™ and HP™ to be A3s, hence examples of 
such manifolds. Combining their revision with Theorem l4.5l one may obtain regularity of optimal 
transportation on any Riemannian submersion of the round sphere simply by showing that no 
minimizing geodesic on the submersed manifold is purely focal. 
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or any of the weaker hypotheses proposed by Loeper [Llj . Notice that the distance 
function is not differentiable anywhere on the cut-locus (at least in CP™ and HP"). 
Otherwise, there is a nonzero gradient Vj,dist(iro, y) of the distance function at y in 
the cut-locus of Xq. Since the domain of exponential is the round ball, moving from 
y along the direction V y dist(a;o, y) should decrease the distance, a contradiction. 
Therefore, the continuity result we obtained also shows that optimal maps stay 
uniformly away from the cut-locus, which enables one to apply the method of Ma, 
Trudinger and Wang [MTWj as done by Loeper |L2j to get higher regularity of 
optimal maps for smoother source and target measures. 

It is natural to expect, but unknown to us, whether every Riemannian submer- 
sion quotient of the round sphere satisfies the strict convexity of the domain of 
exponential map which is the only additional condition required for the method of 
the preceding paragraph to apply (at least for continuity of optimal maps). 

For completeness we mention that for the case of a covering map (Riemannian 
submersion with discrete fibers) of the round sphere, it is known that lifting of 
the measures on B to the total space M can be applied to show the regularity of 
optimal maps using established regularity results [L2j on the round sphere. This 
was discovered independently by Delanoe and Ge in Appendix C of DG . An 
alternative approach (in the same spirit to our discussion above) to this covering 
case has also been given by Figalli and Rifford |FRj . 

6. Sphere is almost positively cross-curved 

In this section we show our final result, namely that the standard round sphere is 
almost positively cross-curved. This represents a significant advance over Loeper's 
discovery that the round sphere satisfies A3s. Its proof will require the following 
elementary lemma. 

Lemma 6.1 (Calculus fact). For < p < ir, the function 

a(p) := sin 2 p + psin p — p 2 (l + cos p) 
satisfies a(p) > 0. Moreover a(p) = if and only if p = 0, n. 
Proof. Reparameterize p :— ir/2 + arcsin(A) by |A| < 1. Then 
a(ir/2 + arcsin(A)) 

= 1 - A 2 + (tt/2 + arcsin(A)) V / l — A 2 — (tt/2 + arcsin(A)) 2 (l - A). 

Define 

a(ir/2 + arcsin(A)) 
HA) := — • 

The assertion holds if 6(A) > 0, for |A| < 1. From 

(1 - X)b'(X) = 2 - A + (| + arcsinA)(2A - 1)/(1 - A 2 ) 1 / 2 

one can check 6(-l) = 0, b'(-l) = 0, and b'(X) > 3 if \ < A < 1. Moreover, 
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for — 1 < A < h, which shows b' (A) increases monotonically in this range. Thus 
b 1 (A) and 6(A) both remain positive throughout |A| < 1, completing the proof. □ 

Theorem 6.2 (Sphere is almost positively cross-curved.). The n- dimensional 
sphere M = S n with the standard round metric (i.e., sectional curvature K = X) is 
almost positively cross-curved (|4.7p . a fortiori non-negatively cross-curved. 



Proof. This theorem follows from the following nontrivial (and tedious) calcula- 
tions. Let's first set up the geometric configuration we are going to analyze. Let 
a; be a point in the round sphere S n of diameter tt. Fix two unit tangent vectors 
q, w G T x S n , \q\, \w\ > 0. For (£R, with \t\ sufficiently small, let r(t) be a line 
in T x S n with f(i) = q, r(t) = 0, \r(t)\ < ir, where /, / denote the time derivatives 
4ff(t), -§£if{t) of a function /. Let x(t) be the c-segment x(t) :— exp x r(t). De- 
note x = x(0), p = \r(t)\ (thus < p < 7r), and f = jT^jj- Let (•,•) denote the 
Riemannian inner product. 
To apply Lemma [231 define 

fai\ u u / dist ('^W) 2 N / \ 

(6.1) H := Hess( ) {w,wj. 

2 x 

. . ,2 

To prove (|4. T|) we will first show —H= —j^H > 0, and then the equality case shall 
be determined. By continuity we may assume < p < ir without loss of generality. 
From a standard Riemannian geometry calculation (for example sec [DL [L2 ), one 
can show that 

(6.2) H = \w\ 2 -IG, 
where 

I := \w\ — (r,w) , G := 1 . 

sin p 

Step 1: reduction to 2-dimensional case. One of the key points of the proof 
is to rearrange the expression of —H in a clever way to enable further analysis. 
Before differentiating H, let's list some preliminary computations in the order of 
complexity. Define a function 

ucosu 

g(u) = : , u G\u,tt[. 



Then, 



where 



<?» - -r—B, g"{u) 
u=p sin p 



P A. 



u=p sin 3 p 



^ 2(sinp — pcosp) p — cospsmp 

p sin p 

Here one can check that A, B > for < p < tt and this will be important later. 
We use these and the identities 

P=(f,q), -(|g| 2 - (r,g) 2 ) 

P 
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to do the following differentiations and rearrangements: 

(6.3) G=-?—A, 

2 sin p 

G=J-B(f,q), 
sin p 

G = -r^-A{f, q) 2 + ^—B{\q\ 2 - (f, q) 2 ), 
sin p psmp 

2 

J = ~(-(r,w)(q,w) + (r,w) 2 (r,q)), 
P 

2 

/= —{A{r,q){r,w){q,w) - 4(r,w) 2 (r, q) 2 - (q,w) 2 + {r,w) 2 \q\ 2 ). 
P 

Key observations here are first, 

G > 

and second, the quantities /, I are independent of the normal component w of w 
to the plane £ C T x S n generated by f and q. Let w\ = w — w 1 be the projection 
of w to E. By separating / = Iw^l 2 + | wi | 2 — (f, w) 2 , one sees 

(6.4) -H = GI + 2GI + GI 

= G\w ± \ 2 -H 1 
> -Hi 

where Hi is the quantity defined by replacing w in (|6.1[) with w\ , thus independent 
of w . Notice that the quantity —H\ becomes identical to —H of r,q,wi viewed 
as tangent vectors of the 2-dimcnsional round sphere that is the exponential image 
of S in the original sphere. This reduces the consideration to two dimensions. 
We will need the following key expression obtained by (|6.3|) and rearrangement: 

(6.5) -m = { [A-£-(f, q) 2 + B(\q\ 2 - (f, q) 2 )] (\ Wl \ 2 - (f, w,) 2 ) 

psmp I sin p 

+ 4(5 - A)((f, q) 2 (f, Wl ) 2 - (f, q) (f, wj) (q, Wi)) 



A{(r, Wl ) 2 \q\ 2 -(q, Wl ) 2 )). 



Note that 



p + p sin p cos p — 2 sin p 

6.6 B - A = — r-!- > for < p < ir, 

p sin p 

as can be checked by taking the fourth-order derivative of the numerator. At 
this point, Loeper's result in |L2] that S n is A3s can be obtained by substituting 
(q, w) = (q, w\) = into the expression (|6.5[) and using the second line of (|6.4[) . 



Step 2: two dimensional case. From (|6.4ll it suffices to show —H\ > 0. From 
now on, we assume without loss of generality the dimension is two, and let f = (0, 1), 
q = (cos6>,sin6>), w\ = (cos ip, sin i/j) in R 2 = T X S 2 , with < 0,if> < 2ir. Let 
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T := tan 6», S := tani/>, -00 <T,S < +00. One checks from (|6T5]) . 

(6.7) -H X = COs2gC ° s2V, p, 

psinp 

where 

o 2 

P = AS 2 - 2(2.8 - A)T S + A-^-T 2 + B - A. 

sin p 

Thus it suffices to show P > 0. P is a convex (since A > 0) quadratic polynomial 
in S 1 with discriminant 

o 2 

D :=4(2B - A) 2 T 2 - AA(A-^T 2 +B- A) 

sin p 

2 

=4(((2S - A) 2 - A 2 -^)T 2 - A(-B - A)\. 
L sin p J 

We show D < (regardless of T), which implies P > 0. Since A(B — A) > 0, D is 
always negative if 

0> (2P - A) 2 - A 2 -^— 



sm p 

= (2B - A + A-?—)(2B -A- A—^—). 

sm p sm p 

The first factor is positive, and the second factor is negative since 

2B-A- A-?— = a(p), 

sm p p sm p 

where 

a(p) :— sin 2 p + psin p — p 2 (l + cos p) 

which is positive from Lemma 16.11 (since < p < tt) . This establishes the desired 
inequality (|4.7|) . 

Step 4: equality case. Let us analyze the cases of equality, to conclude the 
almost positive cross-curvature (|4.7p of S n . We only need to show for < p < tt 
that ~H = holds if and only if the three vectors q, w, f at T x S n are all parallel. 
The necessity is easy to verify. For sufficiency, suppose — H = 0. From (|6.4p . 
w 1 - =0. Thus it reduces to two dimensional case as in Step 2. Now, from (|6.7[) 
and P > 0, cos 6* cos tp = 0. Thus either q or w is parallel to f. In cither case 
examining with (|6.5p shows the other vector is also parallel to f. This establishes 
almost positivity of the cross-curvature of S n . □ 
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